BOUNDEDNESS OF SMOOTH BILINEAR SQUARE FUNCTIONS AND 
APPLICATIONS TO SOME BILINEAR PSEUDO-DIFFERENTIAL 
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Abstract. This paper is devoted to the proof of boundedness of bilinear smooth square func- 
tions. Moreover, we deduce boundedness of some bilinear pseudo-differential operators associ- 
ated with symbols belonging to a subclass of BSqq. 
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1. Introduction 

Let us begin with some important results from the linear theory of singular multiplier opera- 
tors on L P (M). The prototype of a singular multiplier operator is the Hilbert transform, which 
is defined as: 

Hf(x):=-p.v. [ f(x-y)^-, feS(R). 

Or equivalently 

where <S(M) denotes the Schwartz space of functions on K. 

It is well known that H maps L P (R) into L P (R) for 1 < p < oo and L X (E) into L 1 '°°(IR). The 
study of such singular integral operators comes under the "Calderon-Zygmund theory" . 

Let u be an interval in R. Denote by \ u the characteristic function of oj and consider the following 
operator 

ir u f(x) := [ lUOHOe 2 ™^. 
Jr 

Remark that the operator ir^ truncates the frequency to the interval co and if oj = [a, b], a < b, 
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it has the following relation with the Hilbert transform. 



TTu = -{M a HM_ a - M b HM. b ), 

where M is the modulation operator given by M a f{x) := e 2max f{x). Hence, using the bound- 

edness of the Hilbert transform one can easily deduce that the operator tt w possesses the same 
boundedness properties as the operator H. Moreover, the operator norm does not depend on 
the interval oo. 

Observe that for p = 2, using Plancherel theorem we can write 



where uj n , n£Z, are disjoint intervals in M. such that their union is whole of M. The above equa- 
tion and the uniform bound for the operators TT UJn on L P (R) motivate the study of boundedness 
properties for the square function / — > ( ^ |vr Wn (/)| 2 )2 when p ^ 2. At this point we would like 



to remark that the Hilbert transform has one preferred point of singularity, whereas these square 
functions have infinitely many points of singularity. Therefore the study of these operators is 
more delicate. The first result in this direction is due to Littlewood and Paley ([28, 29, 30]). 
They proved that 



Theorem 1.1 ([28]-[30]). Let uj n = (-2 n+1 , -2 n ] U [2 n ,2 n+1 ) 
exist constants c p and C p such that for all f £ we have 



n G 7L. For 1 < p < oo, there 



(1) 



LP( 



< 



<C r 



LP{ 



In literature such square functions are referred as non-smooth Littewood-Paley square functions. 
A little about the proof of inequality (1): First one proves the boundedness of a suitable smooth 
Littlewood-Paley square function. Then vector valued arguments permit to get the right side 
inequality in (1) and the left side inequality is deduced using duality. 

Later, in the year 1967, Carleson [13] considered the non-smooth Littlewood-Paley square 
function associated with the sequence u) n = [n,n + 1], n£Z, and proved the following. 



Theorem 1.2 ([13]). For 2 < p < oo, there exists a constant C p such that for all f G S( 
have 



aw 



(2) 



^K[n,n+1](/)| 2 )' 



<C r , 



LP I 



The smooth analogue of Carleson's Littlewood-Paley square function can be defined by taking a 



smooth function (j) with suppcj) C [0,1] and 4> n (£,) = 4>(£,— n ), ^ £ Z, as / ( ^ |(e 



1-Kin. , 



This smooth operator satisfies the similar LP estimates, i.e., for 2 < p < oo, there exists a 
constant C p such that for all / € 5(K), we have 



(3) 



(>> 2 ™^)*/|^ 



n 



< a 



LP( 



Lp( 



Note that in both these Littlewood-Paley inequalities (inequality (1) and (2)), the sequences of 
intervals have specific properties. In the first case intervals are dilates of each other by a power 
of 2 whereas in the second one they are integer translates of each other. The question for other 
sequences of intervals remained open for quite some time. Finally, in the year 1985 Rubio de 
Francia [36] provided a positive answer towards it in its greatest generality. His result is: 
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Theorem 1.3 ([36]). Let u n be an arbitrary sequence of disjoint intervals in 
2 < p < oo, there exists a constant C p such that for all f E <S(K), we have 



Then for 



(4) 



(£k%(/)l ! 



< C r 



The proof of Theorem (1.3)is quite intricate. The first step towards it is the reduction to the 
case of well-distributed collection and then to invoke some vector-valued Calderon-Zygmund 
theory. 

We now describe some of the existing results concerning bilinear square functions. 

The study of bilinear multiplier operators (or multilinear operators) attracted a great deal of 
attention after the breakthrough of Lacey and Thiele [25, 27] on Calderon's conjecture about 
the boundedness of the bilinear Hilbert transform. The bilinear Hilbert transform is defined for 
f,geS(R) as: 

H(f,g)(x) :=p.v.- / f(x + y)g(x-y) — . 

Or equivalently 

H(f,g)(x) -ill mmV>gn^- V )e 27Ttx ^ ) dCd V . 

Jl Jr 

Lacey and Thiele proved the following. 

Theorem 1.4 ([25, 27]). For the exponents Pi,P2,Ps satisfying 1 < pi,P2 < oo and — + — = 
-r- < I, there is a constant C > such that 

(5) \\H(f, g)\\u3(R) < C||/||iPi(R)||ff||i^2(R). 

We would like to remark that the bilinear Hilbert transform H has the property that 

H(MJ,M a g)(x) = M 2a H(f,g)(x). 

This property is called the modulation invariance or modulation symmetry. A very precise time- 
frequency analysis is required to resolve this symmetry in order to obtain LP estimates for such 
operators. We refer the reader to the works of Gilbert and Nahmod [15, 16] and of Muscalu, Tao, 
and Thiele [34, 35] for the study of operators closely related to the bilinear Hilbert transform. 

Consider the operator 

(6) *M,9){x)= I I K09(v)lUt;-v)e 2 ™^dZd m f,geS(R). 

Jr Jr 

Note that here the one dimensional interval to gives rise to a strip {(£, rj) E M? : £ — n E u} in 
M 2 . So the singularity is along the lines £ — r\ = a and £ — r\ = 6, where a, b are the endpoints of 
the interval us. Using the same arguments as in the linear case, one can conclude the same LP 
estimates for the operator 7r w as for the bilinear Hilbert transform, with bounds independent of 
the interval ui. 

We consider the bilinear Littlewood-Paley square function (/, g) — > ( ^ \^u> n (/, <?)| 2 ) 5 associated 

with a sequence of intervals {u! n }nez- At this point we would like to remark that unlike the 
linear case there is no passage available from smooth square functions to non-smooth square 
functions in bilinear setting. So, in the bilinear case these are two different problems. The 
first result in this direction is due to Lacey [22] . He considered the bilinear analogue of smooth 
Carleson's Littlewood-Paley square function (3) and proved the following: 
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Theorem 1.5 ([22]). Let % £ C°°(R) with supp X C [0, 1]. Define X [n,n+i](0 = nG 
Then, for 2 < p\,p% < oo satisfying — + — = \, there is a constant C > such that 



(7) 



n,n+l] 



The proof of this theorem specifically depends on the exponent 2. Mohanty and Shrivastava [31] 
extended Theorem 1.5 for other exponents p% ^ 2. Their result is: 

Theorem 1.6 ([31]). Let cf) G 5(E). Define <j> n (£) = <j>{£ - n), n G Z. T/ien /or 2 < p ± ,p 2 < oo 
and | < P3 < oo satisfying — + — = i, i/iere is a constant C > suc/i i/iai 



P3 



(8) 



(Ei T ^(/^)i 2 )^ 



< C||/||lpi(k)||3||lp2( 



LP3( 



In the proof of this theorem, the authors bounded the square function by the bilinear Hardy- 
Littlewood maximal operator (studied by Lacey in [23]). That is why the range of exponents 
depends on the one of this bilinear maximal function. By using a particular case of the result 
which we have proved in this paper, we can conclude Theorem (1.6) for the remaining exponents 
P3 G (1, 2] (See § 6). The bilinear analogue of smooth dyadic square function has been addressed 
by Bernicot [8] and Diestel [14]. They proved that 

Theorem 1.7 ([8, 14]). Let ip G <S(M) be non-negative, have support in ^ < |£| < 4 and be equal 
to 1 on 1 < |£| < 2. Define 'ip n (0 '■= ^(2~ n £)> n G Z. Then for exponents Pi,P2,P3 satisfying 
1 < pi,P2 < oo and ^ + ^ = ^ < |, t/iere is a constant C > suc/i i/iai /or a// /, <7 G 5(R) 



(9) 



(^|T^(/, 5 )|^ 



nS2 



< c 



Z,pi( 



LP3( 



Here the authors have linearized the square function using Radamacher functions and showed 
that the resulting operator falls under the setting of Gilbert and Nahmod [15] and then they 
can deduce the desired conclusion. 

The LP estimates for non-smooth bilinear Littlewood-Paley square functions seem to be very 
hard to prove. There is only one result known in this direction, which is due to Bernicot [8]. He 
proved the following. 



Theorem 1.8 ([8]). Let oj n 

b n -i - a n -i and a n+ \ - b n - 
satisfying — + — = — there is a constant C 

V pi P2 P3 

S(M) we have 



- [a n ,b n ] be a sequence of disjoint intervals in R with b n — a n = 
a n — b n -\ for all n G Z. Then for exponents 2 < Pi,P2,Ps < oo 
C{p\,p2,Pz) such that for all functions f,g G 



(10) 



(Ek(/,s)f 



n&L 



< c||/IIlpi(r)II9|Ilp2( 



LP3( 



One can easily observe that the bilinear analogue of non-smooth Carleson's Littlewood-Paley 
theorem (2) is a particular case of Theorem (1.8). The proof of this result is based on time- 
frequency techniques. The author has introduced a new notion of "vector trees, vector size, and 
vector energy" (see §3 of [8] for precise definitions) in his work and proved appropriate estimates 
for these quantities to obtain the desired result. 

The proof of the previous theorem relies on standard time-frequency analysis (used for the 
bilinear Hilbert transforms) mixed with £ 2 -valued arguments to deal with the square function. 
But it is based on a geometric assumption : the intervals co G are exactly well-distributed, 
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they have the same length and are equi-distant. 

We can expect that the result remains true for more general collections O but the analysis seems 
to be very hard. That is why, in this paper we aim to study the smooth square functions, which 
should be easier to study. 

In this paper we have obtained LP estimates for smooth bilinear square functions. Some of our 
arguments are general and can be extended to more general situations. But, we could complete 
the proof of our main result only under the following assumption: 



(11) 



inf | a; | ~ sup |w|. 



Theorem 1.9. Let £l : = (a;) wg ^ be a well- distributed collection of intervals satisfying (11). 
Then for exponents p\,p2,p'^ € [2, oo] satisfying 

1 1 1 

< — = — + — , 

P3 Pi P2 

there exists a constant C, independent of the collection Vl, such that for all f,g G <?(R), we have 

1/2 

< c||/IIlpi(r)IIs , IIlp2(ir). 

LP3 (R) 



\uiefi 



By the counterexample in [8], we know that pi,P2 > 2 is necessary to get boundedness in above 
theorem. 

Let us point out that in this work, our approach allows to consider the case of a collection of 
intervals with equivalent lengths. Whereas in [8], the result relies on a tricky approach based 
on the fact that the intervals are well-distributed (with equal lengths and equi-distant). 

Moreover, we would like to describe another motivation to our work, concerning bilinear pseudo- 
differential operators. We refer the reader to [7, 6] for works of the first author and of Benyi, 
Torres and Nahmod [3, 9]. There, several classes of bilinear pseudo-differential symbols are 
studied (see Section 5 for more details). We are specially interested in the "exotic" class BSq 
which is the following one: a function a E C°°(]R 3 ) belongs to BSqq if 

^|^,e,r?)| <C aAc , 

for all indices a, b and c. This class corresponds to the bilinear version of the linear class of 
symbols S® . For such a symbol, we look for boundedness of the associated operator 



T a {f,g){x):= [ e ix ^f(09(vMx,^ V )d^dr,. 
Jm. 2 



In [12], Calderon and Vaillancourt proved that the linear operators associated with symbols in 
<Sq are bounded on L 2 (IR). A natural question arises for the bilinear operators ? 

In this framework, a first negative answer has been given by Benyi and Torres in [1]: for 
every exponents p,q,r £ [l,oo) satisfying Holder rule, there exist symbols a £ BSq (which are 
x-independent) such that the bilinear operator T a is not bounded from L P (R) x L q (R) to L r (M). 
From this point, a question arises : which extra assumption on a symbol a G -B5qo yields the 
boundedness of the corresponding operator ? 

This work has also the motivation to give some answers to this question for the local-L 2 case. 
We refer the reader to Section 5 for more details concerning existing results. 
The class BSq is invariant under rotation in the two frequency variables. Unfortunately, the 
study of the non-smooth square functions in Theorem 1.8 is based on the time- frequency analysis 
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used for the bilinear Hilbert transforms. It is now well-known that for such operators, there are 
three "forbidden" lines 

{(M, £ = 0}, {(£,??), r? = 0} and {(£,77), £ + r? = 0}. 

So, Theorem 1.8 can be extended to the case where symbols l u (£, — rj) are replaced by symbols 
— tr]) with i ^ {0,-1, 00}. Since we are interested in the class BSq (where these three 
degenerate lines do not play a role), we do not want to use Theorem 1.8. Indeed our new ap- 
proach for Theorem 1.9 will allow us to consider these "degenerate" lines as well as the other 
ones. 

These considerations are also another motivation for Theorem 1.9. We will define classes 
Wq' s (BSq 0) of bilinear symbols (included in BSqq) (see Definition 5.2) and obtain the following 
result for the associated bilinear operators. 

Theorem 1.10. For every 9 G 8 1 and s G (1,2], consider a symbol a G Wq' s (BSq ). Then the 
associated operator T a is bounded in the local L 2 -case, i.e., T a is bounded from L P (IR) x L 9 (IR) 
into L r (R) for exponents p,q,r' G [2, 00] satisfying 

1 1 1 

< - = - + -. 

r p q 

Moreover, the estimates are uniform with respect to 9 G § . 



The current paper is organized as follows. 

2. Notation 

For all m G K and p G (1, 00), we set W m ' p (M) for the fractional Sobolev space on R, defined 
as the set of distributions / G S'(R) such that J m {f) G L P (R), where J m := (Id - A) m/2 . 
For a function / G S(R), the Fourier transform of / at the point £ is defined by 

/(£)= / /(•'•)' '•'"' </•'•• 



With this definition, the inverse Fourier transform is given by / v (£) = (2tt) 1 f(— £). 
For a bounded symbol a, the bilinear operator 



is well-defined and gives a bounded function for each pair of functions /, g in S(K). Moreover for 
for bounded symbols a, the operator T a maps S(K) x S(R) into S'(R) continuously. This justifies 
many limiting arguments and computations that we will perform without further comment. 
The formal transposes, T* 1 and T* 2 of an operator T : S(R) x S(R) — > S'(R) are defined by 

(T*\h,g),f)=(T(f,g),h) = (T* 2 (f,h),g), 

where (■, •) is the usual pairing between distributions and test functions. 

For an interval oo of the real line, we write Xuj f° r a "smooth scaled version" of the characteristic 
function 1^. It means that \u> is a smooth function with support contained in u and for every 
integer i G N, the following estimate holds 



D (i) 



X 



< CM 



To 0, := (w)cjgn a collection of intervals, we associate the following bilinear square operator : 

% :=(/,<?) ^ (yjT x M,g)\ 2 ) ■ 



SMOOTH BILINEAR SQUARE FUNCTIONS 

This can be considered as a smooth version of the more singular bilinear operator 

1/2 



Definition 2.1. A collection := (u) of intervals is said to be well-distributed if 

(12) £ 12* < I- 

wen 

Part 1. Boundedness of bilinear square functions 

This part is devoted to the proof of Theorem 1.9. First we deal with the limiting case when 
r = 2. 

Proposition 2.2. For p,q € [2,oo] satisfying 

1 1 1 

2 P q 

there exists a constant C such that for every well-distributed collection £1 satisfying (11), we 
have 



\\Sa(f,g)\\ L 2 {R) < C||/|| iP(E) || 5 | 



We do not write the proof and refer the reader to Theorem 1.1 in [22]. Indeed in [22], M. Lacey 
dealt with the particular case where f2 := ([n,n + l]) n ez and Xu> are obtained by translations. 
However, we point out that he obtained an estimate of 



^X[«,n+l] (/' 9) 



(uniformly with respect to the symbol employed for X[n,n+i])> see inequality (3.6) in [22]. Then, 
he summed up all these estimates for n € Z and obtained the desired inequalities. So in the proof 
of M. Lacey, we can use different symbols X[n,n+i]- Proposition 2.2 for the case f2 := ([n, n+l]) ne z 
is proved in this way. Then we extend it to a general collection of intervals by remarking that if 
Xj n is a smooth scaled version of lj n with J n C [n, n + 1] and \J n \ ~ 1, then it is also a smooth 
scaled version of lr nin+ i]. Hence, we obtain Proposition 2.2 for any such collection of intervals. 
Since exponents satisfy Holder relation, using the invariance under dilations, we conclude to the 
whole Proposition 2.2. 

Moreover, using invariance under dilation we can replace Assumption (11) by the following one 

(13) — < inf \uj\ < sup \cj\ < 10, 

10 wen u( zq 

where the scale is fixed. 

As a consequence of this, in order to prove Theorem 1.9 it suffices to prove the following : 

Theorem 2.3. Forp,q,r' G (2, 00) satisfying 

1 _ 1 1 

r p q ' 

there exists a constant C such that for every well- distributed collection £1 satisfying (13), we 
have 

\\Ssi(f,9)\\Lr(s.) - C||/IUp(K)lblU?(M)- 
Remark 2.4. Let us point out two improvements. 
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• The proof is based on a time-frequency decomposition where the functions f, g and the 
symbols will be decomposed with wavelets. So we can replace the symbol Xlu{v~€) (defin- 
ing the bilinear multiplier operator n^) by any symbol m^iy], £) supported in {rj — £ £ to} 
and such that 

for all multi-index a. 

• Moreover, we let the line 77 — £ = play a role in the singularity of the symbols considered. 
Indeed we can consider other lines. We refer the reader to Remark 4-14 f or the fact that 
we can consider the lines 

{(£,»?), A 1 £ + A 2 r ? = 0} 

as soon as X\ ^ \i- 

The next two sections are devoted to the proof of Theorem 2.3. 

3. Reduction to the study of combinatorial model sums 
First, we "regularize" the collection Q in the following sense : 
Lemma 3.1. To prove Theorem 2.3, we can assume that the collection Q satisfies 
(14) E 1™ < 1 

for a large enough parameter k > 2. 

Proof. Let S] be a well-distributed collection of intervals and k > 2 be fixed. For each uj £ Q, 
we will consider a "finite Whitney covering", built as follows. For co = [c(cj) — v>c(w) + ^r], 
we define for i = —1,..,N — 1 

e \ M 1 ir * i + 2 i 

uji := c{u) Ml^r, ~jyH- 

So we have a covering oj C Ui=— 1 We associate a partition of the unity consisting of smooth 
functions Xu,i supported in oj{ such that 
JV-l 

• Yl Xu,i = 1 on w 
i=-i 

• For all i and for all n, we have 

The last implicit constant depends on N, but TV" will be later fixed by k. 
So we have the following decomposition : 

N-i 

\ T xM,g)\ < Yl \ T x»x«,Af'9)\ ■ 

i=-l 

The symbol XuXu,i is a smooth cutoff relative to the subinterval Wj. So it remains us to check 
that for all i G {—1, ...,N — 1}, the collection f2j := (co^^^n satisfies (14) for some integer N. 
Then, we estimate the initial operator Sn by the sum of N + 1 operators associated with 
collections Qi verifying (14). Since N will be finite and chosen with respect to k, we can conclude 
the proof of the Lemma. 

We fix an index i € {—1, ...,N — 1} and study the collection fij. We set N > 4k. First, note 
that we have = jt\oj\ so 

2k, 1, 

\K0Ji\ = — M < -\(jJ\. 

1 n ]\T 2 



< c a 

00 



< \uji\- n < \u 
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In addition, 

c(Wi) = c(w) - — + G w- 

Consequently, Kxjj C 2w. Since the collection f2 is well-distributed, we conclude that the collec- 
tion fij satisfies (14). □ 

From now on, we only consider a collection f2 satisfying (13) and (14) for some large enough k 
(which is sufficient due to the previous lemma). 

Since our bilinear operators have modulation invariance property, to prove Theorem 2.3 we 
require the "standard" time-frequency analysis used for the study of bilinear operators, such as 
the bilinear Hilbert transform. We have to decompose both in the frequency and in the physical 
space with the notions of "tiles" and "tri-tiles" . 

For each interval w G O, we use the duality and introduce a trilinear form as follows : for all 
functions f,g,h € S(R) 

(15) (T x M,g),h)= [ /(fr)£(6)£(£3)x w (6-£ 1 R- 

•^l+&+f3=0 

We recall the notion of tiles and tri-tiles (see for example [34, 35]) and adapt the definition to 
our setting. 

Definition 3.2. A tile is a rectangle (i.e. a product of two intervals) I x to of area one with 
1/10 < |7| < 10. A tri-tile s is a rectangle s = I s x to s , which contains three tiles s, = I Si x u Si 
for i = 1, 2, 3 such that 

(16) V^G {1,2,3}, I Si =I s , 



(17) € to Sl +io S2 +uj S3 , 

and such that there is one ( and only one ) interval well satisfying 

(18) V(&,&) eco Sl xoo S2 , 6-&ew 

and \I s \\u\ — 1- Let us enumerate the collection Q := (w n ) n >o- Then for all n > and Q a 
collection of tri-tiles; we will denote by Q n the sub- collection of tri-tiles verifying (18) for the 
interval to = to n . 

The condition (17) describes the fact that the cube to Sl xw S2 xw S3 meets the plane £3), £i + 

£2 + ^3 = 0} (over which the integral in (15) is computed). If a cube does not intersect this plane, 
then the corresponding tri-tile does not play a role and so need not be considered. Moreover, 
we emphasize that as the bilinear symbol is supposed to be smooth "at the scale 1" (since all 
the intervals to € Q have a length equivalent to 1 due to Assumption (13)), we need to deal with 
only those tiles whose space and frequency intervals have lengths equivalent to 1. 
We can cover the time-frequency space 

»x (\J{(Z,vU-veu}) 

where the square function Sq, plays a role by tri-tiles. 
We recall the concept of grid and collection of tri-tiles : 

Definition 3.3. A collection T := {I}j^x of real intervals is called a grid if for all k 6Z 

(19) Yl 11 ^ Co1k ' 

iei 

2 fc-l<|/|<2fe+l 

where Cq is a large enough numerical constant independent of k and the collection I. 
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The classical dyadic grid is a grid in the previous sense. Moreover, a grid has a similar 
structure as the dyadic grid, i.e. at each scale 2 fc , the collection of the intervals (whose the 
length is equivalent to 2 k ) of the grid is a bounded covering of R. 

We remember the notion of collection of tri-tiles (see for example Definition 2.8 in [23]) 

Definition 3.4. Let Q be a set of tri-tiles. It is called a collection of tri-tiles if 

• {I s , s £ Q} is a grid, 

• J := {uj s , s £ Q} U (Ui=l Wsi, s £ Q}) is a grid, 

• for all oj' € J and s £ Q 

3i £ {1, 2, 3}, u> s% C J £ J => V j E {1, 2, 3}, w flj . C a/. 

The last point corresponds to a "sparseness" assumption about the collection J (including all 
the frequency intervals of the tri-tiles), We refer to Definition 2.1 of [33] for more details about 
the notion of "sparseness". 

Now we define the wave packet for a tile. 

Definition 3.5. Let P = I x uj be a tile. A wave packet associated with the tile P is a smooth 
function $p satisfying 

• Fourier support of <5p is contained in ^jW. 

• For all indices i £ N and all M > 0, following estimate holds 

-M 



d« [ e **")-$ P ] ( x )| < [/[- 1 / 2 - 4 (i + |x 



where c(I) denotes the center of the interval I and the implicit constant depends on 
exponents M. 

So, $p is a normalized function in I? (R) , concentrated in space around I and its spectrum is 
exactly contained in uj. 

Remark 3.6. As an instance of the Heisenberg uncertainty principle, the normalization \I\\u\ ~ 
1 for a tile P = I x co is chosen. This ensures that wave packets on P exist. We know that we 
cannot choose a function $p which is perfectly localized in P, i.e. localized in both the time and 
the frequency space simultaneously (which is the case in the Walsh model, see [34] ). However 
we just require some decay for the function $p around the physical interval I, but it raises some 
technical difficulties (see [35]). 

Since in Theorem 2.3 we are only considering exponents strictly bigger than 1, we can use duality. 
Hence in order to prove Theorem 2.3 it suffices to obtain bounds for the following trilinear form 

(20) (f,g,h) -+W T XaJn (f,g)(x)h n (x)dx 

for / £ L p (R), g £ L q (R) and h = (h n ) n a sequence belonging to L r ' (R, Z 2 (N)). 

Then the first reduction is to pass from the "continuous" operators T Xuin to discrete variants 
involving sums of inner products with the wave packets (introduced in Definition 3.5). This step 
is standard and appears essential in order for the phase plane combinatorics to work correctly. 
This reduction could be performed by decomposing functions /, g, h n with wavelets (in terms of 
Fourier series on some intervals scaled with respect to the strips associated with the frequency 
intervals ui £ 0). We refer the reader to [10] and [11] for a precise explanation of such reduction. 
Consequently Theorem 2.3 can be reduced to the following one about model sum operators. 
Indeed, following the ideas of [10] and [11], the trilinear form (20) is equal to an absolutely 
convergent sum of disretized trilinear forms Aq (see (21)) with different collections of tri-tiles 
Q and collections of wave packets. Theorem 3.7 below gives a uniform bound for such model 
trilinear forms and hence similar estimates hold for (20). 
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Theorem 3.7. Let Q be a collection of tri-tiles. For functions f,g € S(R) and a sequence 
h := (h n ) n of smooth functions, consider the following trilinear form : 

(21) AqCf,^)-^ Yl \Isr 1/2 \(f^ Sl )(9^s 2 )(h n ^ sa )\. 

n>0 s£Q„ 

Then for 2 < p, q, r' < oo satisfying 

1 _ 1 1 

r p q : 

there is a constant C (depending only on exponents and implicit constants appearing in Defini- 
tions 3.3 and 3.5) such that for every collection of tri-tiles Q, 

|Aq(/,5,/j)| < C||/|| LP ( R )||5|| L9(R ) 



We further reduce our problem using real interpolation. Since the range for exponents is an 
open set, the strong type estimates (corresponding to Theorem 3.7) are implied by weak type 
estimates using multilinear Marcinkiewicz interpolation. 

Definition 3.8. For a measurable subset E ofM, we write : 

F(E) := {/ e L°°(R), Vx e R, \f(x)\ < l E (x)} . 

Let pi,p2,P3 be positive exponents. We say that a trilinear form A is of weak type (pi,P2,P3) "if 
there exists a constant C such that for all measurable sets Ei,E2,Es of finite measure and for 
all functions f G F(E\), g G F(E2) and sequences h := (h n ) n with YlneZ I 1 2 £ F(E^) we have 

3 

(22) lAC/^/OI^CnW 7 * 

The best constant in (22) is called the weak type bound for the trilinear form A and will be 
denoted by C(A). 

With the help of Marcinkiewicz real interpolation theory applied to the bilinear square functions 
(see the work of L. Grafakos and T. Tao [19] or Exercise 1.4.17 of [17] for bilinear interpolation 
and the one of L. Grafakos and N. Kalton [18] for extension to sub-bilinear operators), weak type 
estimates for exponents varying in an open range yields strong type estimates. Thus, Theorem 
3.7 is reduced to the following one (see Remark 8 of [8] for similar arguments) : 

Theorem 3.9. Let 2 < Pi,P2,Ps < oo satisfy 

1 1 1 

— + — + — = 1. 

Pi P2 P3 

Then the trilinear form Aq is of weak type (pi,P2,P3) uniformly with respect to any finite col- 
lection of tri-tiles Q and associated collection of wave-packets. 

4. Study of these combinatorial model sums. 

To prove Theorem 3.9, we organize the collection Q into sub-collections called vectorized tri- 
tiles and then we study orthogonality properties between them. We emphasize that the following 
estimates do not depend on the collection Q. 




Let us recall the notion of sparseness: 
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Definition 4.1 (Definition 4.4 in [35]). A collection of intervals I := is said to be sparse, if 
land I' are any two intervals in I with \I'\ < \I\ and 10 5 In 10 5 /' ^ <fi, then either \I'\ < 10 9 \I\ 
or I = I' . By extension, we say that a collection of tri-tiles Q is sparse if the collections 

{I s , s£Q} and {u) Sl ,u S2 ,u) S3 , s G Q} 

are sparse. 



We leave it to the reader to check that there exists a finite number K such that every collection 
of tri-tiles can be split into K sparse sub-collections. So without loss of generality, we can assume 
that the considered collection is sparse. 

As in [8], we will require the notion of vectorized tri-tile: 

Definition 4.2. Let s C Q be a tri-tile and j G {1, 2}, we define s J to be the following collection 
of tri-tiles: 

s- 7 := {s' G Q, Sj = s'j] . 

Then, we define for {j, 1} = {1,2} 

s .— U tg ^j i . 

Remark that if t £ s •? then It = I s . 

Definition 4.3. Let j G {1,2} and D := (sj')j be a collection of I -vectorized tri-tiles (with 
I 7^ j)- We say that D is strongly j-disjoint if 

(1) for all i ^ i' and for all s G sf, s' G sf/ 1 , we have Sj Pi s'j = 

(2) and for i / i' and s G sr, s' G i^, if 2u Sj f] 2uj s i / then I s ,, n I Si = 0. 

We now define the quantities size and energy. 

Definition 4.4. Let P be a collection of tri-tiles, f a function, j G {1,2} and I ^ j. We define 

^(/,p) ^supi^- 1 / 2 [ v i(/,$ fl 

For a sequence of functions it is defined as: 

Definition 4.5. Let P be a collection of tri-tiles and h := (/i n )nez be a sequence of functions. 
We define for {j,l} = {1,2} 

\ 1/2 

sl?4\h,P) := supl/.r 1 / 2 ( ^ ]T l^.fc. 

Definition 4.6. Lei P be a finite collection of tri-tiles, j G {1, 2}, / 7^ j, and f be a function. 
We define 

1/2 



eners&(/,P) := sup sup 2 fc [ V |J S | I , 

k&DCP J 



where we take the supremum over all the collections D C P of j-disjoint l-vectorised tri-tiles 
s l G D such that 

2 2k \I s \< |(/.^.)| 2 <2 2fc+2 |I s |. 
We define the vectorized energy for a sequence of functions as follows: 
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Definition 4.7. Let P be a collection of tri-tiles and h := (h n ) n be a sequence of functions. 
Then 

1/2 



energ^h, P) := sup I ^ ^ |(/i n , $ Sj )| 2 j 



where we take the supremum over all the collections D of different tri-tiles. 
In the next section we prove estimates for all these new quantities. 

4.1. Estimates of the quantities "size" and "energy". Let us first obtain required esti- 
mates for the quantities associated with single function. 

Theorem 4.8. Let P be a finite collection of tri-tiles, f be a function and j € {1, 2} and I ^ j, 
then 



(23) siz^(f,P)<\\ JU «, { 

Remark 4.9. The proof given below, can be slightly simplified. Indeed, keeping the Assumption 
(11) in mind, we decompose the frequency plane into cubes of length equal to one. So the 
geometry of such a decomposition is simple. 

The proof, we give here is general and shows how the arguments are related to the geometry of 
the strip. With our approach we are able to use the linear square function associated with the 
collection £1. This proof also holds for an arbitrary well- distributed collection of intervals. 

Proof. As the collection P is finite, we may fix a tri-tile s such that 

1/2 

/ 1 . — - i 2 \ 

(24) 



5«5(/,P) = [4, £ |(/>^;.)) 

\ |s| s 'eV / 



We know that ~s^ 1 is a collection of tri-tiles having the following properties : for each integer n, 
s 1 n P n is a collection of tri-tiles s' whose the frequency part lj s ' x uj s ^ x oj s ^ satisfies 

• |o; s /| = \oj s > 2 \ = \w s '\ ~ \ui n \ 

• for all (£i,£ 2 ) G w s / x uv 2 , £2 - £1 G w n . 

Consequently, due to the property of the collection of tri-tiles (see Definition 3.4) and sparseness, 
we get that 

tt(^nP„) e {0,1}. 

Let denote uj s ^ n the (eventual) frequency cube corresponding to this tri-tile. So we have the 
following description of s l 

~f l = 1J/ S X U s>n . 
n 

Consequently, we get 

1/2 



y S| n>0 



sizeJ(/,P) < I ifrEk/'^x^),) 



Since the functions $>i sX (u s n )j are normalized in L 2 (IR), have Fourier support contained in (oJ s ,n)j 
and have very fast decay around I s , we deduce that 

m>0 1 s ' 

where integer M can be taken as large as we want. We recall that ir^ a n y is the linear Fourier 

multiplier associated with the symbol l( Usn ).. 

Moreover, it is easy to see that there exists a point £0 such that 

(25) (u s ,n)j C Co ± c"LU n 
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for some numerical constant c". 

Indeed, let us fix j = 1 and / = 2 for example. By definition, we know that 

• (w s )2 = (w s „)2 where s n is the tri-tile belonging to ~~{? 2 n P n 

• and if no is the integer such that s belongs to P no; we have 

|(w 3n )i| ^ \uJ n \ ^ |w no | ~ 

where we have used Assumption (11). 

Let ?7o be any point of (uj s )2, then for all n > 0, we have 

(w, B )i C (w s Ji - (w fl „) 2 + (w s J 2 
C -W n + (w s J 2 
C -w n + T] - c(u n ) + cu} n 

Cl] - CUJ n , 

with numerical constants c, d . The other situations (the case of integers j ^ I) can be studied 
exactly the same way. So (25) is proved. 
Hence we conclude that 



si4(/,p)<E 



??? 



-M 



1 



m>0 



1/2 



. n>0 



L 2 (ml s 



The point £o acts as modulation and does not depend on n, so we get 



s~riej(/,P)< 



m 



-M 



1 



m>0 



Im/J 1 / 2 



1/2 



El 7r ± C "-n(^°"/)P 



,n>0 



L 2 (ml 3 ) 



Thus the linear square function relative to the collection (±c"u n ) n >o comes into the picture as 
remarked earlier. Since the collection (w n ) n >o := is assumed to satisfy (14) for some large 
enough k (see Lemma 3.1), the collection (±c"u; n ) n >o is still a bounded covering of R. So we 
can apply Rubio de Prancia's inequality (see [36]) to the linear square function 



1/2 



(26) 



,n>0 



Indeed we will just apply the easy part for p = 2. At this point we would like to point out 
that the appearance of this linear square function is the main reason for considering only the 
"local-L 2 range" for the exponents in the Theorem 2.3. 
Let us now decompose / as follows: 



/ - /l4I s + ^ / 1 2 i '+ 1 / s \2 i, / s 



b>2 

For the first term, with a large enough integer M, we get that 

1/2 



£ 

m>0 



■/?? 



-AI+l 



1 



m 



LW 2 



< 



J~n/P llj ll iP ( 4/ s ) 



< 



L 2 {ml a ) 



For the terms f^-2 b + 1 i s \2 b i s with 6 > 2, we use the same reasoning by noticing the following fact 

22bM 1 



(A 



2 b + 1 I a \2 b I a i ®s 



-26Mi 



-/ i 2 6+i/ s \ 2 ''/ s5 U h — nr^J 
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We remark that the function ) 2A/l & Sj is always a wave packet for the tile Sj (according to 

Definition 3.5). As (1 + %^) ~ 2 b for x in 2 fe+1 / s \2 b / s , the same arguments hold as previously 
and we get an extra factor 2~ 2bMl with a power M\ as large as we want. Consequently we obtain 



that 



size;.(/,p) < ii/iil- ( r) + E 2 ~ 26Ml ii^ii^ W S II 



b>2 



which concludes the proof of (23). □ 

Theorem 4.10. Let P be a finite collection of tri-tiles. For j G {1, 2}, I ^ j and a function f , 
we have (independent of the collection P) 

^^•(/,P)<ll/lb(R)- 

Here we follow the arguments given in [35] and provide the detail of the proof for easy reference. 
We first prove some lemmas. 

Lemma 4.11. Let P be a collection of tri-tiles and set a Sj := (f,<& Sj ). Then there exists 
a collection U of strongly j-disjoint vectorized tri-tiles (st')i and complex numbers c s . for all 
s G Uj Ht 1 such that 



energ 



g(/,p) 



EE 



a Sj c Sj 



and such that for all i 



E 



sesi 



I 2 < 



Proof. Since the collection P is finite, take k and U be an optimizer of the vectorized quantity 
"energy". For all s G Ujsf', define 



Then, 



c Sj :=2- k (Y,\^\r" a s 



EEss = (££k/)2- fe (£|/ s j)-s 



< 2 fc (^|/ Ss |)i 



= energy^/, P). 

The other side of inequality is proved using the second condition. Also the inequality 



£ 



is immediate from the definition of c s . . 



Ic I 2 < - 



□ 



Lemma 4.12. Let U k a collection of strongly disjoint vectorized tri-tiles (sf')i in P and for 
each s G G U, Zei c Sj . 6e a complex number such that 

(27) £ |c s /<A|/ s J, 

/or some A > 0. TTien we /iai>e 



(28) 



£ £ c -A 

i s€sf ! 



< (^£1^1)^ 



L 2 ( 
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Proof. Taking the square on both sides of (28) it is enough to prove that 

E E Kc sfj m Sj ^)\<A^\i Si \. 

Since the Fourier transforms of & Sj and $ s ' are supported in oo Sj and oo s i respectively, the 
term , <Ev.) vanishes if u s . n uj s i = 0. So we need to consider only those terms for which 
u) Sj n oj s i, 0. As the sum in s, and Sj/ is symmetric, we can assume that \I S , \ < \I Si \, which 
implies |u; s .| < |c*v I. 

Using the decay property of $ s , we have 



i 



l^,^)l<^fl + diSt(/ -^ V 



'~|/.|*V w 
< Ijvj^ A dist( j a< ,/ v ; 



. ,- ( 1 + 

1^1 2 



-100 



i *i i 



-100 



since \I S .,\ ~ |J Si | ~ 1. 

Substituting this we only need to prove that 

-100 

(29) 



id' V 1 sJ 7 s ei« s 'e^ t 



l"S^ |<|W / | 



Since all the tiles have equivalent size, we know that \oo Sj \ ~ \uj s i\. Then we estimate using the 
following 



< \c s A 2 + \cA 2 - 



It is enough to estimate the contribution of the first term |c s | 2 , as the second one is similar. 



For a fixed oj s . , the above condition and the sparseness imply that necessarily u s r = u) s . . Hence 

by the definition of j-disjoint vectorized tri-tiles (because s' G sf/ 1 ^ s^ 1 , s',s G would imply 
s = s' due to the sparseness and the fact that j ^ I and the rank 2 of the collection of tri-tiles) , 
the summation over i' and s' have disjoint spatial intervals. As a result, this sum just contributes 
a numerical constant and hence we can estimate the right hand side of (29) by 

E E ki 2 

* self' 

which satisfies the desired inequality because of (27). □ 

Proof of Theorem 4.10: The proof follows applying Lemma 4.11, as we can find a collection of 
strongly j-disjoint vectorized tri-tiles U := {sj'jj, and complex coefficients c s . for all s 6 Ujsf' 
such that 



energy^/, P) 

and 



(/■EEs*. 



E 



,2 < \±Si\ 
s i I ~ I r I 
Z^i \ 1 s i \ 



for all indices i. Applying Cauchy- Schwartz inequality and Lemma 4.12, we get the desired 
result. □ 

Now we prove similar results for the "Z 2 -valued quantities" . 
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Theorem 4.13. Let P be a finite collection of tri-tiles, h := (h n ) be a sequence of functions 
and {j, 1} = {1,2}, then 

(30) siz4 l (h,P) < ||||Mj3(N)|| LooW > 

where the implicit constants do not depend on the functions and the collection of tri-tiles. 
Proof. Let us choose a tri-tile sq such that 




Fix n > and to make computation easy, assume that j = 2 and 1 = 1. We analyze the 
collection U := (w Sj 3) sg ^2,i nPn and prove that it consists of disjoint intervals. Consider two 
different intervals oj S) s and u) s i 3 in this collection. Then by definition, there exist r, t' € 5o 2 
such that (jj 3 \ = u; Tj i and uv 1 = c*V i- Since the collection P n is of rank one, we know that 
r and r' should necessarily belong to different strips : indeed if r, t' belong to the same strip, 
then due to the sparseness and the rank one of the collection "s*" 2 , we deduce that r = r' and so 
uj S) i = ojg',1, which with (uj s ,2 — w S) i)n (uj s ',2 — w s ',i) 7^ (since s, s' belong to the same strip) and 
sparseness yields that s = s'. Hence u St s = uj s > : 3, but it is not possible. However these two strips 
have equivalent widths (since Assumption (13)). Due to the well-distributed property (14) and 
w Ti 2 = u) T ' t 2 (by definition of the 2- vectorized tri-tile), we know that 

d(u) Ti i,(jj T i t x) > d(io p ,ujpi) > k >> 1 

where p denotes the strip of r and p' the strip of r' and k >> 1 as introduced in Lemma 3.1. 
Finally, we deduce that 

(31) d(uj Sj i,uj s ' t \) > k » 1 

Moreover, we know that 

G uj s .i + 0Js2 + oj s ;z and uj s ^ - u s ,l G uj p 

hence 

-W S)3 C U s> i +UJ S A +U} p . 

Similarly, 

Thanks to (31) and \oj p \ ~ |cj p /| ~ 1 << k, the two previous embeddings yield 

d(w S) 3,uv j3 ) > 1 ~ |w S)3 | ~ |w S)3 |. 

Consequently, we obtain that the collection C7 := (^>s,3) S £^j,inP n ^ s composed of disjoint intervals 
of equivalent length. So the corresponding linear square function 

\ 1/2 

E 

seio J ' i nP^ 

is bounded in L 2 (M), uniformly in n. Then we repeat the proof of Theorem 4.8, concerning the 
estimate of the "size" quantity, using L 2 -boundedness of this new square function instead of the 
one given by (26) as before. The same arguments still hold and permit us to achieve the desired 
result. □ 

Remark 4.14. Our aim in this remark is to describe how the arguments depend on the singular 
line (see Remark 2.4). If we replace the symbol Xw( r ? — (defining the bilinear multiplier operator 
Ww) by an y symbol mu(ri,£) supported in {Ai£ — X2I] £ oj} and 



<C a 
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for all multi-index a. Then, we leave it to the reader to check that all the previous results still 
hold (uniformly in X±,X2): indeed they are based on the decomposition of the frequency plane 
into cubes of length equivalent to 1 and the proofs do not use the geometry of the strips (see 
Remark 4-9). 

Only the previous proof uses the structure of the strips and the arguments used degenerate when 
Ai = A2. That is why Theorem 2.3 can be extended to such symbols as soon as X\ 7^ A2. 



We now obtain the following bounds for the quantity "energy" : 

Theorem 4.15. Let P be a collection of tri-tiles and h = (/i n )„>o be a sequence of functions, 
then we have 



(32) 



energ 



k^p)<||i^ 



n|U 2 (N)|| L 2( 



Proof. Let us fix a collection of tri-tiles D. Since each collection D n P n is of rank one, by 
sparseness the collections D n P„ can be assumed to be the collections of disjoint tri-tiles of 
equivalent size. Then it can be easily checked that 

E \( h n,$ Sj )\ 2 < ||MU 2 (R) 
seP„nfl 

which permits us to conclude the proof. □ 



4.2. Proof of Theorem 3.9. We follow the "standard" reasoning employed for this kind of 
time- frequency analysis. So we first obtain the estimate for one vectorized tri-tile and then use 
the combinatorial algorithm to obtain the final estimate for the entire collection of tri-tiles. 

Proposition 4.16 (Tri-tile estimate). Let P be a collection of tri-tiles. Then for each tri-tile 
sq € P, index {j,l} := {1,2}, there exists an implicit constant (independent of the tri-tile sq 
and the collection P) such that for all functions fx, fa 6 S(R) and all sequences fa = (fa, n )n, we 
have 

\,- A.I\.h.fa) :=E E \Is\- 1/2 \(fx,* sl )(fa,*s 2 )(fa,n,$s 3 )\ 

< \I s \size\{fi,P)size\{fi,P)sizef{fa,P). 
Proof. For example, let us assume that j = 1 and I = 2. By definition 

A5+i, 2 (/l,/ 2 ,/3):=E £ |I s r 1/2 |(/l,* S1 )(/2,* S2 )(/3,n,^3)l 
n>0 seiS 1 ' 2 nQ„ 

£ \(fa,*s 2 )\J2 E \(h^s0(hn,*4 



-1/2 



seso 1 



n>0 s'SQn 
s 2 = s 2 



By Cauchy-Schwartz inequality, it follows that 

^^(h,fa,fa) 



( 



1/2 



<\I S0 \- 1/2 E lto.*«>l 

sGso 1 



1/2 



E \(h^s[) 

( 



1/2 



E E Vhn^s 



n>0 s'gq n 
s' 2 =s 2 

\ 1/2 



< ( E ifo.*- 



sup rs 

sGso 1 



t r 1 / 2 



E <A^-i) 



s'sP 
\ s 2= s 2 



E E \(hn,®s 

n>o s 'e4i,2nQ n 
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Using the definition of size, we get 



1/2 



E K/2,^ 2 )| 2 ) <|/^| 1/2 5iS(/ 2 ,P) 



and 



Moreover, we have 



1/2 



E E |(/3,»»*« 
,n>0 s 'eiS 1 ' 2 nQ„ 



< IZ S0 I 1 / 2 ^S^ 2 (/ 3 ,I > )- 



sup |/ s |~ 1/2 

sGio 1 



f ^ V2 

E I (A.*.* ' 



s'eP 



<size?(/i,P). 



This completes the proof. 



□ 



We now prove a combinatorial algorithm which will allow us to organize the whole collection 
of tri-tiles into sub collections in such a way that we have control over the different quantities 
associated with these sub collections so that we can sum them up to get the desired result. 

Proposition 4.17. Let j E {1,2} be fixed and I ^ j. Let P be a collection of tri-tiles and f a 
function such that for some integer d E Z 



(33) 

Then we can decompose P 
(34) 

and the collection P 2 
(35) 



«*4(/,P) <2- d i^g.(/,P). 
pi y p2 SQ ^ e co ii ec n on p 1 satisfies 

4(/. pl ) < 2- d " 1 i^^(/,P) 
~Si l, i)i consisting of I -vectorized tri-tiles Sj'^ , is j- 
U <2 2d . 



■disjoint and satisfies 



E 

i 



Proof. We follow ideas of Proposition 12.2 in [35]. Let us denote the energy E := energy^- (/, P) 
and consider the case I = 2 and j = 1 (the other cases can be similarly treated). We initialize 
with a collection D to be the empty collection. 
We consider the set of all tri-tiles s C P satisfying 



(36) 



E 

s'e 



-in 



> 4" 



2~ d E 



II, 



If there are no tri-tiles obeying the previous condition, we terminate the algorithm. Otherwise, 
we select s± among all such tri-tiles. We add tri-tiles it''- 7 into the collection D and remove 
them from the collection P. Then we repeat the algorithm until there is no tri-tile left satisfying 
the selection criteria. This terminates in finitely many steps as the whole collection P is finite. 
In the process, we have constructed a collection (sj)j of tri-tiles. We set P 2 := Ujst Zj and 
P 1 :=P\P 2 . 



T>2 



IS 



We claim that the selected tri-tiles (sj)j satisfy : the collection of Z-vectorized tri-tiles s, 
strongly j-disjoint. It is clear from the construction that for i < i', tri-tiles s £ Sj' and s' E s^ 1 
are different. In fact we know even more that Sj ^ s'- (else s' would have been removed with 
Sj'' J at the ith step). By sparseness we deduce that Sj n s'j = 0. Now suppose, on the contrary, 
that we had tri-tiles s 6 Sj' and s' E sf/ 1 for i ^ i! such that 2co s j n 2u s ij ^ and L s i C L Si = I s . 
By symmetry, we can assume that i < i! . Since — the sparseness assumption would 

imply that u). 



8 ,3 



which is not possible since s' would also be removed with and 
cannot be found by the algorithm at the step i! > i. We thus arrive at a contradiction, which 



s ;3 ' 
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proves the strong j-disjointness of the vectorized tri-tiles. 

It remains us to check property (35) (since (34) is obvious by construction), which corresponds 
to 

i 

Since the /"-vectorized tri-tiles Si 1 are proved to be j-disjoint, by using the definition of "energy" 
(see Definition 4.6) with (36) and (33), it follows that 

2- M £ 2 ]T|/ Si |<£ 2 , 

i 

which proves the desired inequality. □ 

Now, it remains to prove a similar algorithm for the quantities associated with the sequence of 
functions. Indeed, this part is far more easy, and is proved in the following proposition. 

Proposition 4.18. Let h := (h n ) n be a sequence of functions, {j,l} = {1,2} and P be a finite 
collection of tri-tiles such that for some integer d G Z 

(37) siz(%\h,P) < 2~ d energlj 3 (h,P). 
Then we can decompose P = P 1 |J P 2 where the collection P 1 satisfies 

(38) size^ l (h,P 1 ) < 2~ d ~ 1 Sner^ 3 (/i, P) 
and P 2 is a collection of vectorized tri-tiles P 2 = (st J '')j such that 

(39) ^|/ Si |<2 M . 

i 

Proof. We follow the previous ideas. Let us denote the energy E := energy^ (/i, P). We initialize 
D a collection of tri-tiles to be the empty collection. 
We consider the set of all tri-tiles sCP satisfying 

(40) Yl E Yk^sO^ >^ l {2~ d E) 2 \I s \. 

n s'€~£j> l nP n 

If there is no such tri-tiles, we terminate the algorithm. Otherwise, we choose a tri-tile si among 
all tri-tiles satisfying (40). We add sp'' to the collection D and we remove it from the collection 
P. Then we repeat the algorithm and we construct a sequences of tri-tiles (sj)j. When this 
algorithm is finished, we have constructed a collection P 2 := UjSj J '' and P 1 := P \ P 2 . 
It remains us to check property (39), which as before, is a direct consequence of the fact that D 
is a collection of different tri-tiles by construction. □ 

We can now complete the proof of our main result. 
Proof of Theorem 3.9: 

The arguments are now routine. We refer the reader to [32, 35] for precise arguments. We will 
just give a sketch of the reasoning. Fix functions /i, /2 and a sequence of functions fa = (/3, n )ru 
then we want to estimate Aq(/i, /2, fs). 

Thanks to Propositions 4.16, 4.17 and 4.18, we can iterate the different combinatorial algorithms 
in order to get a partition of the initial collection Q 

Q = U Q d , 

where for each d S Z and for all {j, /} = {1, 2} we have 

si^J(/,-,Q d ) < 2-fegyj(/ i ,Q), 

and 

siz4 I (/ 3 ,Q d ) < 2- d eHe^ 3 (/ 3 ,Q). 
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In addition, the collection Q d can be covered by a collection of vectorized tri-tiles for some 
collection such that 

(41) £|/ s J<2 2d . 

ieD d 

Proposition 4.16 yields 

A Q rf(A,/2,/3)< II min{2- d eliirg^^(/ j ,Q),s^(/ j ,Q)} 

min 1 2~ d {energy 3 (/ 3 , Q) } , size^' 2 (/ 3 , Q) , size 2,1 (/ 3 , Q) j 
< 2 2d J] min{2- d elieTg^^(/,,Q),s"izt^(/ J ,Q)} 



b',0={i,2} 



mm 



2~ d {elieT#3(/3, Q)}, s4^' 2 (/ 3 , Q), ^e 2 3 >\f 3 , Q)} . 



Then we can compute the sum over d and get the desired inequality : 
Aq(/i,/ 2 ,/ 3 )<EV(/i»/2 1 /3) 



W}={h2} 

eherg^ 3 (/ 3 , Q) 1 ~ 9a (s^3' 2 (/ 3 , Q) + s^ 2,1 (/ 3) Q)f a , 



for every exponents 9i G (0, 1) with #i + #2 + #3 = 1 (see for example Proposition 4.3 in [34]. 
Now we invoke Theorems 4.8, 4.10, 4.13 and 4.15 with the choice 0% = 1 — ^ to conclude that 
for all fi 6 F{Ei) we have 

A Q (/l,/ 2 ,/3) < |£l| 1/pi |£ 2 r /P2 |£3| 1/P3 , 

which is the desired weak- type (pi,p 2 >P3) estimate for Aq and it is independent of the collection 
Q. Hence we conclude the proof of Theorem 3.9. □ 



Part 2. Applications 

5. Applications to bilinear pseudo-differential operators 

Let us first recall classes of bilinear pseudo-differential symbols. Two main types of x- 
dependent classes of symbols have been studied in the literature. One is the Coifman-Meyer 
type class BS™ S (M), < 5 < p < 1, m € K, of symbols satisfying estimates of the form 

(42) \<%dp c v a(x,Z, V )\ < C aM (l + |e| + \ v \)™+ s "-P( b +°\ 
for all indices a, b, c. 

The other type corresponds to classes which are denoted by BS™ S . g(M), < <5 < p < 1, mgl, 
— 7r/2 < 9 < 7r/2, and consist of symbols satisfying 

(43) \%%%,v(x,Z,t,)\ < C aAc;6 (l + \v -t a nm) m+8a - p{b+c) 

(where for 9 = tt/2 the estimates are interpreted to decay in terms of 1 + |£| only). Both the 
classes can be seen as bilinear analogs of the classical Hormander classes S™ S (M.) . 

It is now well-understood that the operators in BSiq are examples of certain singular integral 
operators and fit within the general multilinear Calderon-Zygmund theory developed by Grafakos 
and Torres [20]; see also the work of Kenig and Stein [21]. So their boundedness properties in 
Lebesgue spaces are well-known. 

The general classes BS™ S . S with x-dependent symbols were first introduced in [3]. For m > 
and (p, S) = (1,0), their boundedness properties were obtained by Bernicot [6, 7]. Also some 
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results were obtained by Bernicot and Torres [9] for the exotic classes BS^.g. These class of 
bilinear pseudodifferential symbols allow to build functional calculus (see [7, 5]). 

As explained earlier in the introduction, this section is devoted to the study of the bilinear 
operators associated with symbols belonging to the exotic class BSq = BS® g. As pointed 
out, this class does not depend on the angle 9. 

Remark 5.1. First we point out that the behavior on the product of modulation spaces (in- 
stead of Lebesgue spaces) of bilinear operators associated to a symbol belonging to BSq is 
well-understood, see [2] . Results on some products of Besov spaces also follow from well-known 
embeddings. 

In addition, we refer the reader to [5] for the following result: every symbol a £ BSq define a 
bilinear operator which is bounded from I? x W s,oc to L 2 for sufficiently large s. 

Considering boundedness from the product of Lebesgue spaces, we know from [1] (Proposition 
1) that extra assumptions on the symbol are necessary (due to some counterexamples). We also 
consider some smaller classes of bilinear symbols and show that the associated bilinear operators 
are bounded. 

This approach was already treated in [1] where the authors add the following assumption: for 
every integer a > 



(44) 

and 
(45) 



sup / \\d?m(x,£, -)| 



G?£ < oo 



sup 



\d%m(x,-,r])\\ L2 di] < oo. 



In [1], following a bilinear version of the Calderon and Vaillancourt reasoning, Benyi and Torres 
proved that a -BSpQ-symbol satisfying these two previous conditions gives a bilinear symbol 
which is bounded from I? x L 2 to L 1 . Moreover, using the Wigner transform, same results are 
obtained in [1] for a symbol a verifying for all i,j G {0, 1} 

sup ||3|^<r(x, ' 



L 2 



These two results hold in the setting of multi-dimensionnal variables. 

In this current work, we will define other assumptions and prove the boundedness in the local 
L 2 -case for the assocoated bilinear operators. Mainly, we assume that our symbol is very smooth 
along the x-variable but we would like to require less regularity in the frequency plane. 

Definition 5.2. Let 9 € S 1 be a unit vector, we define 9 1 - to be the orthogonal vector in M 2 . 
We say that a symbol m : M 3 — > K belongs to the class w}' s (BSq ) if m £ C QC (M 3 ) satisfies : 
• the symbol m G BSo.o, for all index a,b,c £ N 



d^m{x^ri) 



< oo 



• the symbol m belongs to W ' s in the direction 9 : for all index a € N 



sup 



sup 



d°m(x, .9 + td 



+ 



d£(V \t >v )m(x, .9 + t9 ± ),9) 



< oo. 



In this case, we set WmW-,,,!^ , uqQ to be the maximum of all the above constants. 

So a symbol belonging to W ' s (BSq ) is a BSq -symbol which satisfies an extra Sobolev regular- 
ity only along the direction 9 in the frequency plane. In [1] (see (44) and (45)), the assumptions 
were weaker in space (since we require less regularity in the variable x) and stronger in frequency 
(regularity in the two frequency variables is required). 
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Theorem 5.3. Let a : R 3 — > R be a BSq^ symbol such that there exist a direction 9 G S 1 
with s/20 £ {(1,-1), (-1,1)} and s G (1,2] such that a G Wg' s (BS$ ). Then the operator T a 
is bounded in the local L 2 -case, i.e. it is bounded from L P (R) x L 9 (R) into L r (R) for every 
exponents p, q, r' G (2, oo) such that 

1 _ 1 1 

r p q 

Remark 5.4. The end-point estimate r = 2 remains true for x-independent symbols a. 



Proof. The case of an x-independent symbol: 

For simplicity, we will assume that IMI^M^go j = 1- We also assume that a depends only on 

the variable A := ((£,r]),0), i.e. c(^,rj) = m(A) for some function m. The idea of the proof is 
the following: we will decompose the symbol m with elementary functions in order be bring the 
previous square functions into the picture. 

Consider a smooth function x supported on [— 1 , 1] such that for all £ G R 

We decompose the symbol as follows: 

pez p&L 

Note that the symbols m p are well-adapted to the interval [p — l,p + 1] with a norm bounded 
by 

rp+i 

sup \m(t)\ < / (|m(A)| + |m'(A)|)dA. 
For each integer k, we denote D k the following set : 

D k := |p, 2 k < J (|m(A)| + \m'(X)\)d\ < 2 k+1 
Since \m\ + \m'\ G £ S (R), we deduce that 



1 : = IHIw^ > || H + KHI J,. = J (|m(A)| + \m'(X)\) s drj 

>E r +1 (KA)| + |m'(A)|rdA 

>E(/7 |m(A)l + K(A)|dA ) 

(46) >^22 ks ($D k ). 
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So only for k satisfying 2 k < 1, the contribution to the above sum is non null. We now estimate 



< 



Ei T -p(/>s)i 



< 



< 



E 

2 fc <l 



E \ T m P {f,g)\ 



L r ( 



/2 



feez 

2 fc <l 



1/2 



E I T ^p(/'5')| 



where we used Cauchy-Schwartz inequality. Then, from what we described previously and 
Theorem 2.3, we know that 



1/2 



E l T "w(/>9)| 



~ II/IIlp(R)II5||z,9(]R)2 



L'( 



As a consequence, we conclude that 



2 k <l 



Finally using (46), we obtain 



\TM 



< 



< 



LP(R)\\9\\Li( 



Lp(R)\\9\\li( 



1/2 



fcez 

\2<=<1 



1/2 



where we have used that s < 2. 



Thus, we have proved the desired result for those x-independent symbols which depend only on 
the above variable A := ((£,77),$). We leave to the reader to verify the general case when the 
symbol a depends on the two variable £, r]. We produce the same reasoning: by decomposing 
(£,77) into the basis (9,9 ) and then work on the main variable ((£,T}),9). All the estimates 
remain true since the other variable ((^,rj),9 ± ) does not play a role and the assumptions are 
uniform on it. 



The case of an x-dependent symbol: We refer the reader to [6] for details concerning the 
"general" principle which allow us to obtain boundedness for rr-dependent symbols by knowing 
the boundedness for corresponding ^-independent symbols. This consists of two steps, we briefly 
recall it. 

For x-independent symbol, the previous decomposition only uses square functions associated 
with intervals of lengths equivalent to one. Using the ideas of [6], we can obtain "off-diagonal 
decay" at the scale 1 for the bilinear operators associated with x-independent symbols. More 
precisely, we can prove that for every interval I of length 1 and m a bilinear multiplier belonging 
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to Wq' s (BSq ), there exists S > (as large as we want) such that 



viL ]T > 



>(/, g)(x)\ r dx 



(47) 




\f(x)\ p dx 



1/p 



g(x)\ q dx 



1/9 



This improvement is obtained in exactly the same way as described in the previous sections with 
the following slightly sharper estimate. We first decompose the tiles s according to the distance 
between the space intervals d(I s ,I). Then, it is important to get fast decay according to this 
quantity. To achieve this, we require new bounds for the "size" quantities, so we need to prove 
(instead of Theorem 4.8) that for j = 1,2 and P a collection of tri-tiles, we have 

-N \ l l 2 



size$(/,P) <sup( 1(1 + 



Is\ JR V Us 



\f{x)\ z dx 



for a large enough integer N. We again leave it to the reader to check that the proof described 
for Theorem 4.8 allows us to get this stronger estimate. Then this improvement gives us the 
off-diagonal decay for the bilinear operator T m (as detailed in [6]). 

Then, by Sobolev inequality, we can extend these local "off-diagonal" estimates for re-dependent 
symbols a belonging to Wq' s (BSqq) (see [6] for details). As a consequence by summing these 
off-diagonal estimates (for / = [n, n + 1], n £ Z), we conclude the boundedness of the associated 
bilinear operators on the whole space. □ 

The main result of this section now can be obtained as a corollary to the previous theorem. We 
recall that the range of exponents in the local-L 2 case is invariant under duality. 

Theorem 5.5. Let m : M 3 — > M be a BSo : q symbol. Suppose that there exist a direction 9 and 
s £ (1,2] such that m £ Wq' s (BSq ). Then the operator T m is bounded in the local L 2 -case: 
from L P (M) x L q (M) into L r (R) for every exponents p,q,r' £ (2,oo) such that 

1 

r 

Moreover, the estimates are uniform with respect to 6 £ S 1 . 



1 1 
- + -• 

p q 



Proof The previous theorem gives us the result for all the directions 9, except one degenerate 
direction. The idea is then to use duality in order to get around this technical problem. 
Let us first remark that for an x-independent symbol m belonging to Wq' s (BSq ), we have the 
following descriptions of the two adjoint operators : Z^ 1 := T m * and := T rn * with 



m{(£,r)) = m(-£ - 77,77) 
ml(£,,rj) = m(f,-J7-£)- 
Hence, m £ Wq' s (BSq ) is equivalent to £ W Q ' s (BSq ) which further is equivalent to 
* 2 eW^ s (BS° 0fi ) with' 



rn 



n*2\ 



cot(0) + cot(6»* 1 ) = -1 and tan(0) + tan(0* 
where cot and tan are naturally defined in § . 

So the off-diagonal decay obtained for multipliers associated with the class Wq' s (BSq ) can be 

transferred to the dual classes W^ s (BSq ) and W^ s (BSq ). Then, as described in the previous 
proof, we can get boundedness on the whole space for the x-dependent symbols belonging to 
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the same classes. 

Since the previous theorem gives the result for all the direction except one, this reasoning relying 
on duality allows us to get the result for this specific direction. Moreover, it is easy to check 
that the estimates are uniform with respect to the directions 6 £ S 1 . □ 

6. Proof of Theorem 1.6 for remaining index 1 < p 3 < 2 

We provide here the reasoning for the remark we made in the introduction concerning the 
validity of Theorem 1.6 for remaining index 1 < p 3 < 2 using our main result Theorem 1.9. 



Proof. Let (j) G 5(R) be as in Theorem 1.6. Consider a smooth function x supported on 
such that for all f G R 

i = 5>(£-p). 



■1,1] 



We decompose the function (ft in a similar fashion as in previous section : 



HO 



\ 



Using Minkowski's inequality we have 

El^n(/^)| 2 ] 



LPS I 



< 



E 




LP3( 



LP3 (R) 



where (j> p>n (0 = 4> P (C ~ n). 

Note that for each p£Z, the symbol 4> p is well-adapted to the interval [p — l,p + 1], hence we 
can apply Theorem (1.9) to the square function associated with 4> Pt n, to conclude 



Yl \ T <i>P,n(f>9) 



< CI 



LP1(M)\\9\\LP2( 



LP3< 



where exponents pi,P2, and p 3 satisfy the Local-L 2 condition. For each p G Z using the decay 
of cf) we know that for N G N, there exists a constant Cn such that 

Cn 



y p||£°°( 



< 



(1 



Hence summing over p, we get the desired result. 
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